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Abstract 
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instability of hadronic matter was found, as a direct prediction of Statistical Bootstrap. 
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1. Introduction 

Quantum Chromodynamics is unquestionably the fundamental theory of strong inter¬ 
action. However, the nonperturbative aspects of QCD belong still to a held of intense 
investigation. In fact, most of the properties of the hadronic world cannot be extracted 
yet from hrst QCD principles and the partition function of interacting hadrons produced 
and thermalized in high-energy collisions can only be determined within specihc models of 
varying degree of limitation. 

In this paper we employ Statistical Bootstrap, including volume corrections for the hnite 
size of hadrons, in order to describe the hadronic phase. The virtue of this description is 
associated with the fact that, among Statistical models of hadrons, it is the only one which 
predicts the onset of a phase transition in strongly interacting matter and therefore it is 
compatible with the properties of the QCD vacuum at high temperatures [1]. The aim 
of this work is to pursue the search for this compatibility with QCD even further, asking 
whether Statistical Bootstrap of hadrons is consistent with the existence of a critical endpoint 
in strongly interacting matter, at high temperature and hnite (nonzero) baryonic chemical 
potential. The existence of such a critical point in the phase diagram is required by QCD in 
extreme conditions [2,3] and it is the remnant of chiral QCD phase transition [2,3]. 

The basic ingredients in our approach are (a) the hadronic partition function extracted 
from the equations of Statistical Bootstrap and (b) the equation of state of the quark-gluon 
phase given by recent QCD studies on the lattice [4-7]. Our principal aim from the matching 
of these two descriptions is to trace the formation of a critical point in the quark-hadron 
phase transition with a mechanism compatible both with Statistical Bootstrap and Lattice 
QCD. 

In [8,9] a similar scheme was pursued but with the quark-gluon partition function med¬ 
itated by the MIT bag model. This model is a crude approximation leading either to an 
ideal gas of quarks and gluons [8] near the critical line or to a gas of weakly interacting 
quarks and gluons with perturbative terms in the sector of two light flavours [9]. The only 
non perturbative effect in this model comes from the pressure of the vacuum through the 
phenomenological bag constant. In order to have an adequate description of QCD matter in 
this approach the use of lattice equation of state, derived from hrst principles is an essential 
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improvement. To this end we have employed, in this work, a realistic partition function of 
QCD matter which has become available recently as a result of a breakthrough achieved in 
lattice QCD studies allowing for computations at nonzero chemical potential [4,7]. As far as 
the hadronic phase is concerned we continue to use in this work the pressure partition func¬ 
tion but with the volume Laplace conjugate variable ^ as an active thermodynamic variable 
and not fixed at the zero value. These calculations, within a three flavour bootstrap scheme, 
are used here for the first time. Despite all these alterations we shall show that the existence 
of a critical point at finite chemical potential persists and we shall try to narrow the area of 
its location. 

In Section 2 the principles and basic equations of Statistical Bootstrap, leading to the 
partition function of hadronic phase, are briefly summarized and the appearance, within 
this framework, of a thermodynamic instability is discussed in detail. This instability is the 
origin of the formation of a critical point in the bootstrap matter. In Section 3 the partition 
function of the quark-gluon phase is extracted from recent lattice calculations of the pressure 
of QCD matter [4]. In Section 4 the above two descriptions of strongly interacting matter are 
exploited in a search for a critical point in the phase diagram, the location of which is fixed by 
a set of equations (29-33). Finally, in Section 5 our conclusions are presented whereas in the 
Appendix certain technical points are discussed concerning the evaluation of the quark-gluon 
partition function (part A) and the complete form that acquires the equation of maximum 
hadronic pressure (part B). 

2. Thermodynamic Instability in Statistical Bootstrap 

The main attribute of the Statistical Bootstrap (SB) [10-13] is that it describes thermo¬ 
dynamically interacting hadronic systems. Generally as far as its thermodynamic behaviour 
is concerned a system of strongly interacting entities may be considered as a collection of 
particles with the explicit form of the complex interaction among themselves. Another way 
to deal with the problem is to assume that the strongly interacting entities form a greater 
entity with certain mass. Then, if this mass is known, the greater entities may be treated 
as non-interacting and the simple thermodynamic description of an ideal system of parti¬ 
cles may be applied. The SB adopts the second approach and assumes that the strongly 
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interacting hadrons form greater clusters, called “fireballs”. The problem then is moved to 
determine the mass of these clusters or to evaluate their mass spectrum f{m?)dim?, which 
is equal to the number of discrete hadronic states in the mass interval {m, m + dm}. The 
solution is based on the assumption that the hreballs may consist of smaller hreballs with 
the same mass spectrum or “input” particles which may not be divided further. The integral 
bootstrap equation then reads [14-17] 


■H^k( 


oo T / n \ 

b,s,...) = 9b,.s,...B{p‘^)6o{p‘^ - ml ) + — / 5^ Ip-^Pi] 

input term 

n \ / n \ n 

B{pl)f{pl, bi, Si, .. .)d^Pi , 


( 1 ) 


{bi} \ i=l / {si} \ i=l / i=l 

where ~ — b{p>^ — ^ ) 9 {pq). Equation (1) also imposes conservation of 

four-momentum p and any kind of existing additive quantum number (baryon number b, 
strangeness s, etc.), through the Kronecker function Sk, between a hreball and its con¬ 
stituent hreballs. The masses mb,s,... are the “input” particle masses which constitute the 
smaller hreballs that can be formed and gb,s,... are degeneracy factors due to the spin of the 
“input” particles. The correct counting of states of a hreball involves, apart from its mass 
spectrum which is of dynamical origin, a kinematical term B{p‘^). This term is related to the 
appropriate volume of a hreball which is considered to be carried with it, = ^p^. The 
imposition of the requirement that the sum of volumes of the constituent hreballs has to be 
equal to the volume of the large hreball, as well as, momentum conservation lead to the fact 
that all hreballs possess the same volume to mass ratio. This ratio can be connected to the 
MIT bag constant, ^ The term B{p^), appearing in eq. (1) is then 

2V^Pn 2Vm 
(27r)^ (27r)^ 

Eq. (2) corresponds to a specihc choice for B. Since in (1) the mass spectrum f is 
accompanied by the term B, the bootstrap equation (apart from the input term) remains 
unchanged if f and B are redehned in such way so 


B{p^ 


Bf = B't 


( 3 ) 


Every such transformation leads to diherent thermodynamic description of the system of 
hreballs, since the relevant quantity for this description is only the mass spectrum f' and not 
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the term B'. These transformations are not nniqnely determined in the general bootstrap 
scheme, allowing for varions versions of the model [15]. 

The bootstrap eqnation can be simplihed after performing a series of Laplace transfor¬ 
mations to acqnire the form 


.fdi, {A}) = 2G(ft {A}) - exp[G(;3, {A})] + 1 . (4) 


In the last eqnation G{j3, {A}) is the Laplace transform of the mass spectrum with the 
accompanying kinematical term 
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and v 9 (/ 5 , {A}) the Laplace transform of the input term 


( 5 ) 


{A}) = Y Y. T E E T' 

b'=—oo q'=—oo s'=—oo | 5'|=0 

( 6 ) 

The masses mbqs\s\ correspond to the masses of the input particles, which in this paper will 
be all the known hadrons with masses up to 2400 MeVs, the gbqs\s\ are degeneracy factors due 
to spin and the A’s are the fugacities of the input particles. Here we have used the extended 
version of SB [17], where the states are characterised by the set of fugacities relevant to 
baryon number Xj,, electric charge Xq, strangeness A* and partial strangeness equilibrium 7 * 
(in the following we shall refer to this set of fugacities as {A}, for short). The last fugacity 7 * 
is related to the number of strange quarks plus strange antiquarks [18], whereas the fugacity 
As is related to the number of strange quarks minus strange antiquarks. The introduction 
of 7 s is necessary for the accurate theoretical prediction of the experimentally measured 
hadronic multiplicities. 

The bootstrap equation dehnes the boundaries of the hadronic phase since it exhibits a 
singularity at the point 

{^}) = ln4-1 . (7) 




From the physical point of view this singularity is connected with the behaviour of the 
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mass spectrum as the mass tends to infinity 


f\m\ {A}) 2C'({A})m-“-^ exp[mp*{{\})] , (8) 

where f3 = T~^ and j3* corresponds to the inverse maximum temperature. After a certain 
point, as temperature rises, it is more preferable for the system to use the given energy in 
producing more hadronic states (since their number rises exponentially) than in increasing 
the kinetic energy of the already existing states. For our specific choice of B in (3) the 
exponent a = 4. 

In order to turn to the thermodynamics it is necessary to consider the fireball states in an 
external volume This volume must be distinguished from the physical volume which is 
carried by each fireball. The partition function of the pointlike interacting hadrons is then 




E E E E 7. 


b'l 


b'=—oo 
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AO A \s'\)e-^^P^dp^ = 
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J (27r)^ • (9) 

Every choice for the mass spectrum in (3) (which leads to a certain exponent a in (8)) leads 
to a different partition function and so to a different physical behaviour of the system. The 
usual SB choice was a = 2, but more advantageous is our choice a = 4. With this choice 
a better physical behaviour is achieved as the system approaches the hadronic boundaries. 
Quantities like pressure, baryon density and energy density, even for point-like particles, no 
longer tend to infinity, as the system tends to the bootstrap singularity. It also allows for the 
bootstrap singularity to be reached in the thermodynamic limit [19], a necessity imposed by 
the Lee-Yang theory. Another point in favour of the choice a = 4 comes from the extension 
of SB to include strangeness [14,15]. The strange chemical potential eqnals zero in the 
quark-gluon phase. With this particular choice of a, ps acquires smaller positive values as 
the hadronic bonndaries are approached. With the choice a = 4 the partition fnnction can 
be written down and for point-like particles it assnmes the form 


lnZ,A(Y-^/3,{A}) 


A-BV 

x^G{x, {\})dx = Vfip, {A}) . 


( 10 ) 


5 



For a = 4 the input term acquires the form 

27r^pB ^ . 

where the index “a” runs to all groups of hadrons, each of which is characterised by the same 
set of fugacities (e.g. Kaons with electric charge Q = 1, N and A Baryons with Q = 0, etc.), 
“i” runs to all hadrons in the same group with different masses and B is the energy density 
of the vacuum (MIT bag constant). 

By including corrections due to the hnite size of hadrons (Van der Waals volume correc¬ 
tions) the repulsive part of the interaction is taken into account. The partition function for 
the real hadron gas is written [19,20] as follows 


‘.;3.{A}) = i: 


2(V®** — p/4B)^p^ 




In the above relation the subscript -|- indicates that each single bracket has to be non¬ 
negative, avoiding, thus, the negative contributions to the volume. The four momentum is 
p^ = Pi and because of its presence the integrations no longer factorize. The factorization 
property can be recovered through the grand canonical pressure partition function [19,20] 


fOO 

7,(13, dVe-f''Z(V,l3,{\}) , 

JO 


(13) 


which is the Laplace transformed partition function with respect to volume. Provided that 
the thermodynamic limit limy^oo(l/V) In Z(/5, V, {A}) exists, the integral in (13) converges 
for the values 

? > ?o(/3. {A}) = Jim [lliiZ(;5,V^,{A})] . (14) 

V—^oo 1 _ y 

If we are constrained to values ^ there is no need to employ Gaussian regularization 
[20] to evaluate (13) for ^ < ^q. Then the pressure partition function acquires the form 


7r(^,/3, {A}) 


1 

e-/(/9 + e/(45),{A}) ’ 


(15) 


where / = In Zpt/V. From (15) it is evident that the value of ^ at the thermodynamic limit 
(14) can be obtained from 

eo = /(/9 + W(4i?),{A}). (16) 
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We have, also, to note that for ^ 7 ^ 0 the critical temperature of the hadronic state at zero 
density, Tqhg, does not depend only on the MIT bag constant, B, as is the case in [14-17], 
but on as well. 

The density and the pressure P of the thermodynamic system can be obtained through 
the pressure grand canonical partition function (15) 


{A}) = A 


df{f3 + ^/m,{\}) 

d\ 


1 - 


1 df{f3+^/m,{\}) 


1 -1 


AB 


df3 


where A is the fugacity corresponding to the particular density, and 


PHG{^,P,{X}) = -f{(3 + ^/{4B),{X}) 


1 - 


1 df{P + ^/{AB),{X}) 


-1 


(17) 


(18) 


4B dl3 

Though volume is no longer an active variable of the system it can be calculated for 
given baryon density and ub (evaluated through (17)) and baryon number < B > which is 

a conserved quantity. The volume would be retrieved through the relation 

<B> 


<V>= 




(19) 


With the use of SB in order to describe interacting hadronic systems we can trace the 
possibility of a phase transition. The study of the pressure-volume isotherm curve is then 
necessary. When this curve is calculated a region of instability is revealed. In fact, this 
curve has a part (near the boundaries of the hadronic domain) where pressure decreases 
while volume decreases also (see Fig. 1). Such a behaviour is a signal of a first order phase 
transition which in turn can be mended with the use of a Maxwell construction. 

This behaviour is due to the formation of bigger and bigger clusters as the system tends 
to its boundaries in the phase diagram. In that way the effective number of particles is 
reduced, resulting, thus, to a decrease of pressure. This is the basic mechanism that will 
produce a first order transition at lower temperatures and a critical point at hnite density. 
To show that this instability in the P — V curve is the result of the attractive part of the 
interaction included in the SB we shall calculate a similar curve using the Ideal Hadron Gas 
(IHG) model with Van der Waals volume corrections (repulsive part of interaction). The 
logarithm of the partition function of IHG (corresponding to (10)) is 


\nZptjHG{V,P,{X}) 

V 

27r^/5 


= VfptiHGi.fi, {A}) = 

Xl[Aa({A}) -h Xfi{X])~^]'^gs.im^iK2ifims.i) , 

a i 


( 20 ) 
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where ga.i are degeneracy factors due to spin and the index a runs to all groups of hadrons 
described by the same set of fugacities. This function can be used in eq. (15) to calculate 
the Ideal Hadron Gas (IHG) pressure partition function in order to include Van der Waals 
volume corrections. The result is that the pressure is always found to increase as volume 
decreases, for constant temperature, exhibiting no region of instability and so no possibility 
of a phase transition. 

The comparison of SB with the IHG (with volume corrections) is displayed in Fig. 1, 
where z/q is the normal nuclear density z/q = 0.14 fm~^. In both cases (SB or IHG) the 
constraints < S >= 0 (zero strangeness) and < B >= 2 < Q > (isospin symmetric system, 
i.e. the net number of u and d quarks are equal) have been imposed. Also strangeness is 
fully equilibrated which accounts to setting = 1. 


3. The partition function of Quark Matter 

Having a description of the hadronic phase at hand, it is necessary to proceed with the 
thermodynamic behaviour of the quark-gluon phase. The QGD equation of state at hnite 
temperature and baryon density calls for non-perturbative methods of approach. Lattice 
calculations have been performed but the power of such methods is limited by the sign 
and the overlap problems [21]. The overlap problem is treated with the reweighting method, 
called the Glasgow algorithm [22]. In [23] the method of imaginary chemical potential is used 
to solve the sign problem. In [24-26] various new methods have been used to tackle with 
the sign and/or the overlap problems. Especially in [26] the overlap problem is eliminated 
completely. In [27] lattice Taylor expansion is used around ps = 0, allowing the exploration 
of the phase transition of QGD at the low density regime. 

Lattice calculations of the pressure of the quark-gluon state have been performed at hnite 
chemical potential in [4,6], using an improved reweighting technique. These publications 
include calculations for rather heavy u, d quark masses. The mass of the u, d quarks is 
about 65 MeV and the strange s quark 135 MeV [4,6]. The calculated pressure of the quark- 
gluon phase at pb = 0 is plotted against the ratio of temperature to the transition 



temperature of quark matter at zero baryon chemical potential T/Tc in Fig. 2 of [4], The 
temperature Tc will be denoted as Tqqgp in the following. The results of this graph are 
extrapolated to the continuum limit by multiplying the raw lattice results with a factor 
Cp = 0.518 [4]. 

The lattice calculations for hnite chemical potential are summarised in Fig. 3 of [4], where 
the difference of pressure at non-zero chemical potential and the pressure at zero chemical 
potential (Ap/T"^ = [P{fi ^ 0,T) — P{fi = 0,T)]/T^) is plotted against T/Tc. Again the 
results of this graph are extrapolated to the continuum limit by multiplying the raw lattice 
results with a factor = 0.446 [4]. In this graph hve curves are plotted which correspond 
to baryon chemical potential of 100, 210, 330, 410, 530 MeV. 

With the use of Figs. 2, 3 in [4], it is possible to calculate in principle the pressure of 
the quark-gluon phase at any temperature and baryon chemical potential. The pressure is 
important, because knowledge of the pressure is equivalent to the knowledge of the partition 
function of the system in the grand canonical ensemble 

^'^Zqgp{V,T,^b) =-;^P{T,iib) ( 21 ) 

In order to have a complete description of the dependence of the pressure on the temper¬ 
ature and the chemical potential we use two sets of htting functions. For constant chemical 
potential the pressure as a function oiT/Tc is htted through 


f\x) 


(l\ 






+ 


02 






( 22 ) 


where ai,bi,Ci,di, fi {i = 1, 2) depend on ps, while for constant temperature the correspond¬ 
ing £t of the pressure as a function of ps is given by 


g{x) = a + 6exp(ca;'^), 


(23) 


where a,b,c,d depend on the temperature ratio T/Tc. The fitting procedure has to be 
performed in a self-consistent way and subsequently it is straightforward to evaluate the 
partition function, as well as its derivatives with respect to hb and T at any given point 
(Ti,/iB i). In particular, to evaluate physical observables connected with the partition func¬ 
tion and drive numerical routines the partial derivatives of the pressure up to second order 
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with respect to temperature and fugacity have to be evaluated. These derivatives are then 
given in part A of the Appendix. 

In Fig. 2 we have reproduced the quark-gluon pressure as a function of the temperature 
for constant baryon chemical potential. The squares are points directly measured from the 
graphs of Fodor et. al and the lines represent the calculation with the hts which has been 
performed on these points, via eq. (22). Fig. 3 is a graph similar with Fig. 2, but we have 
focused on the area which is useful for our calculations, the area where the matching with 
the hadronic phase will be performed. Fig. 4 is a reproduction of the Fodor et. al quark- 
gluon pressure as a function of the baryon chemical potential for constant temperature. The 
necessary hts have been performed with the use of eq. (23). 


4. The critical point in the phase diagram 

After developing the necessary tools to handle the thermodynamic description for the 
quark-gluon phase as it is produced from the lattice, we can search for the possibility of a 
quark-hadron phase transition. 

First we have to deal with the free parameters that exist in our models. In [4,6] where 
the results of the pressure is presented the transition temperature of the quark state at zero 
density Tq (called from now on Tqqgp) is not hxed. In [28], however, the QCD critical 
point is studied with quark mass input values closer to the physical ones and a zero-density 
temperature Tqqgp = 164 ± 3 MeV. Therefore in what follows, in order to limit the free 
parameters existing in our scheme, we choose 

Tqqgp = 164 MeV . (24) 

As far as the hadronic phase is concerned, an upper bound for the parameter Tq hg can 
be hxed at the value 183 MeV. This temperature allows for the best matching of the strange 
chemical potential between the hadronic and the QGP phase [15]. So we shall set 

Tq hg < 183 MeV. (25) 

The fact that Tq hg and Tq qqp acquire different values does not imply a contradiction. At 
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/iB = 0 the strongly interacting system belongs to the crossover regime where the quark and 
hadron phases are indistinguishable. Therefore, the zero-density parameters (Tq qgp,To hg) 
do not correspond physically to a distinct transition between the two phases; they simply de- 
hne the zero-density intercept of the extrapolated critical line to small values of the chemical 
potential, beyond the critical endpoint of the system at nonzero density. As a result, in the 
region = 0, 164 MeV< T < 183 MeV) of the phase diagram, defined by the boundary 
values (24) and (25) a sharp thermodynamic separation between the bootstrap and QCD 
phase may not exist. Although the physics of the system in the crossover regime is not yet 
fully understood, it is a plausible assumption (implied by (24), (25)) that the two phases are 
distinguishable only in the domains: /is = 0, T 164 MeV (hadrons) and T 3> 183 MeV 
(quarks, gluons). 

Turning now our attention to we adopt (14) in order to have always a real pressure 
partition function. For simplicity we also choose to have ^ = const, for every set of (T, {A}). 
This means that (14) has to be valid for every set of thermodynamic variables. Since the 
value of ^ at the thermodynamic limit, ^q, depends on the choice of these variables we 
have to locate the specific set that gives us the highest value of .^o- For this reason we 
calculate on isotherms for different values of fis (fixing the remaining chemical potentials 
in order to fulfil the constraints < B >= 2 < Q > and < S >= 0). It is found that for 
constant temperature .^o rises as a function of fis- Then we calculate .^o on the maximum 
value of hb allowed for each temperature, i.e. on the critical curve (fixing accordingly the 
rest of the chemical potentials). It is found that on the bootstrap critical line .^o rises as 
the temperature is increased. Thus the greatest value of corresponds to T = Tq hg and 
consequently {A} = {1}. So in order to have a real pressure partition function for a constant 
value of ^ all over the space of our thermodynamic variables it suffices to require 

e>eo(ToHG, {A} = {!}). (26) 

Finally as a consistency requirement on the thermodynamics of lattice QCD [28,29] and 
bootstrap matter we impose the constraint 

Tcr.p. < Tq qgp ■ (27) 

Then, if the values for the free parameters are chosen, within the above constraints, one 
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may calculate for a specific temperature the pressure isotherms of Hadron Gas and QGP. 
Assuming that the baryon number is a conserved quantity to both phases, the equality 
of volumes is equivalent to the equality of baryon densities and so the connection of the 
isotherms of the two phases is possible through the relation 

< VhG > = < Vqgp - — - = - HG = QGP (28) 

HG ^B QGP 

The graph of the pressure-volume isotherm can be drawn using the plot of pressure against 
the inverse baryon density (see eq. (19)). Then at the point where the isotherms of the two 
phases meet we have equal volumes for equal pressures. 

Tracing the point where the isotherms of two phases meet, we hnd that at a low temper¬ 
ature the intersection of QGP and SB pressure-volume isotherms takes place at a location 
where the Hadron Gas pressure is decreasing while volume decreases. The resulting pressure- 
volume curve includes an unstable part which has to be repaired through a suitable Maxwell 
construction. This curve includes a region where a hrst-order transition takes place. 

As the temperature increases, there exists a value for which the QGP and SB isotherms 
meet at a point where the Hadron Gas pressure has a maximum. In that case no Maxwell 
construction is needed and since this point is located at hnite volume or not zero baryon 
density (equivalently not zero chemical potential) it can be associated with the QGD critical 
point. As temperature rises more, the resulting pressure-volume isotherm also increases while 
volume decreases without the need of a Maxwell construction and the situation belongs to 
the crossover region. 

A graph that summarises the situations met in the pressure volume isotherms of hadronic 
and quark systems in the neighbourhood of the critical point is Fig. 5. In this hgure the 
hadronic isotherms have been calculated for parameters that fulhl the constraints (25)-(27) 
(To hg = 172 MeV, MIT bag constant ~ 136 MeV, ^ = 177.1 ■ 10"^ MeV^) while the 
quark-gluon isotherms are related to Tqqgp = 164 MeV. For the lower temperature isotherm 
a Maxwell construction is needed to remove the instability from the resulting curve. The 
horizontal line dehnes a partition with two equal surfaces (shaded) and represents the hnal 
pressure-volume curve after the completion of the Maxwell construction. At the temperature 
T = 162.06 MeV the two curves meet at the point of maximum hadronic pressure and so 
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a critical point is formed at finite volume. At a greater temperature the plotted pressure- 
volume isotherm corresponds to the crossover. 

To locate the critical point numerically with the use of the lattice partition function, for 
given parameters Tq hg and Tq qgP) the conditions have to be determined for which the SB 
pressure is equal to the QGP pressure at the same volume, corresponding to the maximum 
SB pressure. Setting the factor of partial strangeness equilibrium 7 * = 1 a hadronic state 
is characterised by the set of thermodynamic variables (T, A„, A^, A*), while a quark-gluon 
state evaluated on the lattice [4] is characterised by the two variables (T, A^). The u and 
d quarks are characterised by the same fugacity — ^'q — To evaluate the 

unknown variables we have to solve the following set of non-linear equations 


5b(T, A„, Ad, As) — Ub QGpir, Ag) 

( 29 ) 

Psb{T-, Ad, As) = Pqgp{T, Ag) 

( 30 ) 

dPsB{Ti, Xu, ^d, ^s) ^ 

dV 

( 31 ) 

( ^ )sB (T; Xu, Ad, As) = 0 

( 32 ) 

(T, Xu, Ad, As) — 2 {Q (T, Xu, Ad, As) = 0 

( 33 ) 


Eq. (29) accounts for the equality of the baryon densities of the two phases which is equivalent 
to the equality of volumes, since the baryon number is a conserved quantity. Eq. (30) is the 
equality of pressures of Hadron Gas and QGP. Eqs. (29), (30) determine the point where 
the two pressure curves meet. Eq. (31) requires that the meeting point of the two phases for 
a certain temperature is equal to the point which maximizes the Hadron Gas pressure and 
so this meeting point is the critical point. The form of this equation is discussed in detail 
in the Appendix. Eq. (32) imposes strangeness neutrality in the hadronic phase. Eq. (33) 
imposes isospin symmetry to the hadronic system, demanded in order to compare with the 
lattice studies. 

The area in the (T, /i^) plane which gives solutions for the critical point compatible with 
the constraints (24)-(27) is depicted in Fig. 6 . The line ^ = ^omax represents solutions for 
the critical point with the requirement that ^ is set to the thermodynamic limit value at 
T = Tqhg- Since the value of ^ at thermodynamic limit is calculated through (16), we first 
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determine ^ and the MIT bag constant B, for given Tq hg, by solving the following two 
equations 

e = /(1/To HG + e/4T, {A} = {1}; T) . (34) 

^{1/To hg + e/4T, {A} = {1}; 5) = ln4 - 1 . (35) 

The last equation determines the critical hadronic curve at zero density. We, then, insert 
the calculated values of ^ and T, for the given value of Tq hG) in fhe system of eqs. (29)-(33), 
which is solved to extract the thermodynamic variables that locate the critical point on the 
phase diagram. All the points on the right of the ^ = /omax curve on the T — hb plane are 
compatible with the requirement (26). 

The compatible area is also limited by the curve of constant Tqhg, according to (25). 
Setting Tqhg and choosing a value of ^ that fulhls (26) we hst solve (35) to determine B. 
Then we solve the system (29)-(33) to draw the line of constant Tqhg- 

In Fig. 6 we have also drawn the line T = 164 MeV, which excludes solutions that do 
not fulhl (27). The shaded area represents all the points that qualify as a solution for the 
critical point, compatible with eqs. (24)-(27). 

Recent lattice QCD studies offer, apart from the quark-gluon pressure which has been 
a basic ingredient in our approach, important results on the existence and location of the 
critical point itself. In [29] the critical point is found to reside at Tcr.p. = 160 ±3.5 MeV and 
Hb = 725 ± 35 MeV, with Tqqgp = 172 ± 3 MeV. These calculations have the drawback 
that have been performed with u and d quark mass which has a value about four times the 
physical value. Improved calculations have been performed in [28], where the light quark 
masses have decreased by a factor of 3 down their physical values. The critical point is 
found now (with Tq qgp = 164 ± 3 MeV) to be at Tcr.p. = 162 ± 2 MeV and fip = 360 ± 40 
MeV, a value which is considerably reduced with respect to the previous one. This point is 
depicted on Fig. 6 with the full star and it falls completely inside the compatible domain 
of the critical point, according to our calculations. In [27] the critical point has also been 
determined using Taylor expansion around = 0 of the 3-flavour QCD equation of state. 
The critical point is found to reside at /i^ = 420 MeV. In [30] a method is exhibited about 
how to locate the critical point using imaginary values of fip- This study also shows the 
sensitivity of the critical point on the strange quark mass. 
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In Figs. 7-8 we illustrate a representative solution for the critical point as well as for 
the QGP-hadron transition line from those which are included in the shaded area of Fig. 6 
on the (T,/is) plane. The chosen critical point is associated with Tq hg = 172 MeV, while 
the rest of the parameters are those of Fig. 5. It is located at Tcr.p. = 162.1 MeV and 
/is = 218.7 MeV. The full circle represents our solution for the critical point. This circle is 
the endpoint of the solid thick line representing the bootstrap calculation of the maximum 
hadronic pressure which is close to the first order critical line. 

In Fig. 7 we also compare our solutions for the critical point with other calculations. The 
full stars represent solutions of lattice QCD. The points depicted as “lattice reweighting I and 
II” are the points from [29] and [28] respectively and the one depicted as “Taylor expansion” 
is taken from [27]. The solid thin line is the first order transition line evaluated in [29] 
with large u, d quark masses. The dotted line represents the crossover in this calculation. 
The two slashed lines represent the magnitude of the slope d^T/d/i^ taken from lattice 
Taylor expansion and which indicate the phase transition line according to [27]. Theoretical 
predictions on the location of the critical point from various models [3,31-36] are depicted 
as open squares in Fig. 7. The labels adopted are those of Table I in [21]. The open circle 
in Fig. 7 is the prediction of the critical point within the statistical bootstrap (.^ = 0) from 
[8], where a simplihed partition function (MIT bag) for the quark phase was used. 

In Figs. 6 and 8 we compare our solutions for the critical point with the freeze-out 
points from different experiments. In Fig. 6 we depict freeze-out points from NA49 at 
158 AGeV [37] for systems of different size {C + C, Si + Si, Pb + Ph) that fall inside 
the compatible domain of the critical point, if the errors in their determination are taken 
into account. These experiments are interesting since they could trace critical fluctuations 
associated with a critical point of second order [44]. In Fig. 8 the comparison is extended to 
a larger number of experiments [38-43]. On the same graph we have also depicted the curve 
oi < E > / < N >= 1 GeV [45] (reproduced from [46]) that hts freeze-out points which 
are spread to a wide region of the phase diagram. It is evident that our calculations set the 
critical point to a location easily accessible by experiments, especially by GERN/SPS. 
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5. Conclusions 

Statistical Bootstrap presents a more accurate description of the hadronic phase than 
the ideal Hadron Gas, since it includes in a self consistent way the attractive part of the 
interaction among hadrons through the mass spectrum, as well as the repulsive part through 
Van der Waals volume corrections. This interaction is crucial to investigate critical phe¬ 
nomena in connection with the state of quark-gluon plasma. Among the predictions of the 
bootstrap model is the limitation of the hadronic phase and the forming of an instability in 
the pressure-volume isotherm near the hadronic boundaries. This instability is connected to 
a hrst order quark-hadron phase transition and the existence of a critical endpoint in the 
strongly interacting matter. 

A more realistic pressure diagram of the quark-gluon phase is available from lattice 
calculations, despite the fact that unphysical values of the light quark masses are still involved 
in these studies. From these results the lattice partition function of the quark state, as well 
as, all the necessary derivatives can be calculated, allowing the evaluation of any physical 
observable. 

The joining of the SB and the lattice partition function for the haronic and the quark 
state respectively, allows for the determination of a critical point at hnite baryon chemical 
potential which can be related to the critical point of QCD. 

More recent lattice calculations [28] drive the position of the critical point to smaller 
values of baryon chemical potential as the values of u, d quark masses approach their physical 
values. It is interesting that the current location is situated in the (T,/i^) plane in a region 
easily accessible by the freeze-out conditions of experiments at the CERN/SPS. 

Setting the free parameters in our model in a way to fulhl certain constraints we are left 
with a compatible domain in the (T, plane for the location of the critical point. Recent 
lattice calculations [28] drive the critical point within the domain of our solutions. 

In a previous work [9,8] a similar solution was found with the use of a simplihed partition 
function for the quark-gluon system, based on the MIT bag model. Therefore, the basic 
mechanism in our approach for the formation of a critical point in the strongly interacting 
matter is not associated with the details of the quark-gluon partition function but mainly 
with the instability of hadronic matter revealed by Statistical Bootstrap. In particular the 
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local maximum of pressure in the P — V diagram (Fig. 5) of hot hadronic matter lies in the 
origin of the formation mechanism of the critical endpoint. 
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Appendix 

A. The htting procedures on curves of constant temperature and chemical potential 
allow us to evaluate the derivatives with respect to T for constant [xb and with respect to 
for constant T. Physical observables, however, are given as derivatives with respect to 
temperatnre for constant fngacity or with respect to for constant T. The evalnation of 
the latter (for the pressure) is given by 
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In eqs. (A.4), (A.5) where the 2nd partial derivative with respect to two different variables 
of P appears, the pressnre is considered as a fnnction of these two variables. 

B . After choosing specihc values for the parameters ^ and B, the reqnirement of eq. (31) 
has to be fnlhlled for a certain temperature T = Ti and in the presence of two constraints 
9i — { ^ )SB (^1’ A(i) Ag) 0 and g 2 = { B (^i) A^, A^) 2 ( Q )sb (^i’ Aci; A^) 
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0. Since the maximum of pressure is found for a certain isotherm, the temperature Ti may 
not be considered in the following as an active variable. The same is true for ^ and B, since 
the maximum pressure is evaluated for constant values of these parameters. With the above 
considerations and eq. (19) we may write 
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Similarly for the constraint g 2 we have 

_ _ 

dXddXu dXgdXu dXu 
Eqs. (B.2) and (B.3) may considered as a system of two equations which can be solved to 
determine dXd/dXu and dXg/dXu 
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Eqs. (B.4) and (B.5) may now be inserted to eq. (B.l) to give 
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Eq. (B.7) is the form of eq. (31) for the maximum hadron pressure in a certain isotherm. 
The main contribution comes from the term SPsb/ dXu, as it is verihed by comparing the 
maximum pressure found in the present work using the full equation (B.7) with the solution 
found using only the first term in [9,8]. 
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Figure Captions 

Fig. 1 Isotherm pressure-volume curve for SB and IHG (both with Van der Waals volume 
corrections using the pressure ensemble with the same value of .^). The SB curve is drawn 
for Tohg = 172 MeV ~ 136 MeV). 

Fig. 2 The pressure of the quark-gluon state divided by T*^ versus the ratio T/Tqqgp, for 
constant baryon chemical potential. The lines from bottom to top correspond to gradually 
increasing values of fip- The squares represent direct measurement from the Figs. 2, 3 in [4] 
which depict the lattice calculation and the lines indicate our hts to these points. 

Fig. 3 Similar graph with Fig. 2. The pressure of the quark-gluon state is divided by Tq gcp. 
The plot focuses on the region where the matching with the hadronic state will take place. 
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Fig. 4 The pressure of the quark-gluon state divided by Tqqqp versus the baryon chemical 
potential, for constant values of the ratio T/Tqqcp. The squares represent direct measure¬ 
ment from Figs. 2, 3 in [4] which depict the lattice calculation. The lines indicate our hts to 
these points. 

Fig. 5 Three isotherm pressure-volume curves for Hadron Gas (using SB) and QGP phase 
(using the lattice pressure of [4]). The low temperature isotherm needs Maxwell construction, 
the middle temperature isotherm develops a critical point and the high temperature isotherm 
corresponds to crossover. 

Fig. 6 Domain in the (T, fip) plane (shaded area) which is compatible with solutions for the 
position of the critical point, after imposing the constraints (24)-(27) on the free parameters. 
There are, also, displayed the freeze-out points from NA49 at 158 AGeV [37] and the lattice 
calculated critical point from [28]. 

Fig. 7 The critical point domain of Fig. 6 and a representative solution of the critical 
point (full circle-“SB/II”) at the (T,/is) plane, as well as the first order part of the quark- 
hadron transition (solid thick line). The stars represent the lattice calculated critical points 
[27,28,29]. The line labelled “lat. rew.” represents the phase boundaries from [29] and the 
lines labelled “Tayl. exp.” represent the phase boundaries from [27]. The open squares 
are theoretical calculations for the critical point [3,31-36]. The open circle (“SB/I”) is the 
calculation of critical point using bootstrap model in [8]. 

Fig. 8 Gomparison of the positions of our critical point solutions presented in Figs. 6-7, with 
the freeze-out points from different experiments [38-43]. It is, also, displayed, with slashed 
line the < E > / < N >= 1 GeV freeze-out curve of [45], reproduced from [46]. 


23 










QGP ' ^ 0 QGP 








Critical Temperature, T (MeV) 

C/X * |J* 



0 100 200 300 400 500 

Critical Baryon Chemical Potential, (MeV) 

X3 C/X * L/* 










Temperature, T (MeV) 



0 200 400 600 800 1000 1200 1400 1600 


Baryon Chemical Potential, (MeV) 





Temperature, T (MeV) 



Baryon Chemical Potential, Pg (MeV) 




